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gradient, buoyant heat flux and maximum mixing depth have 
been considered. Model predicted pollutant concentrations 
have been compared to airborne sampled observations. The 
effect of coastal turbulence not observed by the single 
point meteorological measurements made onboard ship greatly 
affects the advection and diffusion of a plume as it moves 
onshore. Additional measurements/predictions particular 
to the coastal area will have to be incorporated into the 
model for it to accurately predict the onshore movement of 
pollutants. 
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ABSTRACT 



The Ris>a National Laboratory , Roskilde, Denmark, atmos- 
pheric puff dispersion model has been tested for an atmos- 
pheric-marine environment. Phis three-dimensional model 
simulates the release of Guassian pollutant puffs and pre- 
dicts their concentration as they are diffused and advected 
downwind by a horizontally homogeneous, time-dependent wind. 
Atmospheric characteristics such as turbulence intensity, 
potential temperature gradient, buoyant heat flux and maxi- 
mum mixing depth have been considered. Model predicted pol- 
lutant concentrations have been compared to airborne sampled 
observations. The effect of coastal turbulence not observed 
by the single point meteorological measurements made onboard 
ship greatly affects the advection and diffusion of a plume 
as it moves onshore. Additional measurements/predictions 
particular to the coastal area will have to be incorporated 
into the model for it to accurately predict the onshore 
movement of pollutants. 
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INTRODOCTION 



The downwind transport and distribution of atmospheric 
pollutants from an isolated source over land or water has 
become an important environmental factor in today’s society. 
The need to understand the distribution of smoke, unpleasant 
or potentially harmful foreign gases and perhaps radioactive 
debris from a nuclear powerplant accident are becoming more 
and more essential for industral operations and construction 
planning. The dispersion of such atmospheric pollutants is 
commonly modeled by a standard Guassian plume model which 
computes one-hour average plume characteristics. 

The Meteorology Section of the Ris0 National Laboratory, 
Roskilde, Denmark, has recently developed a puff model for 
prediction and simulation of atmospheric pollutant diffu- 
sion. The model considers individual puffs of pollutants 
with specific release rates that are advected by a horizon- 
tally homogeneous wind over a grid. The wind input may be 
either the measured wind from a single point, a spatial 
average or a wind simulation. The model simulates the 
instantaneous plume characteristics by adding a group of 
puffs, growing in size, as they advect with the wind. A 
Guassian plume model, on the other hand, provides a time 
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averaged concentration pattern based on a single time aver- 
age wind vector. In the puff model, the plume advects with 
a time series of actual wind data. Thus, the puff model is 
able to predict time varying concentration distributions in 
actual changing wind conditions, making it an appropriate 
tool for dynamical computations of downwind dispersions of 
pollutants. 

A basic comparison of a puff model simulation and a typ- 
ical plume is illustrated in Fig. 1. Looking from above, 
the instantaneous behavior of a plume being advected from a 
source by the wind is shown. The outer cone-shaped contours 
represent the outer limit of the plume boundary and are 
identical in both Figs. 1 (A) and (B) . 

Fig. 1(A) shows an instantaneous depiction of an actual 
plume. The long-term average plume concentration is shown 
on the extreme right as a smooth curve with a maximum on the 
central axis. Also shown is the instantaneous plume concen- 
tration considered realistic but is of such a short time 
scale that it cannot be predicted or easily measured. 

The puff model prediction is depicted in Pig. 1 (B) . The 
circles show the boundaries of individual puffs of pollu- 
tants released from the source. These puffs are advected 
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PLUME BEHAVIOR 
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Figure 1. Instantaneous Behavior of a rypical Plume and a 
Series of Puffs from a Puff Model 

and diffused downwind by a frequently updated wind. The 
long term average concentration prediction of the puff model 
is expected to be identical to the long term concentration 
of Fig. 1 (A) . The short term average pollutant prediction, 
a Gaussian curve shown on the extreme right, is not com- 
pletely realistic but is a reasonable approximation to the 
instantaneous plume concentration profile. 
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The purpose of this thesis is to evaluate and test the 
general characteristics and capabilities of the Ris0 puff 
model. The adjustable input data will be varied and pre- 
dicted results for various input combinations compared. A 
preliminary comparison will also be made between predicted 
results and data collected from a coastal region using 
observed meteorological forcing data as model input. 
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II, RI^ sail model 



A. GENERAL CHARACTERISTICS 

The Ris0 Puff Model is a three-dimensional computer 
model used for the prediction and/or simulation of the dif- 
fusion and advection of atmospheric pollutants. The puff 
model technique is to simulate a plume with Guassian shaped 
puffs with specified release rates within a specified grid. 
The initial size of the puffs is normally one meter in diam- 
eter although this can be easily adjusted. The amount of 
material in a puff is the release rate times the elapsed 
time between puffs. Therefore, a long elapsed time between 
puff releases results in a higher initial puff pollutant 
concentration than a short time interval. This should not 
normally be of concern if an adequate balance is maintained 
between grid size, advection speed and puff release rate. 

The location of the puffs on the grid is determined by 
computing their movement for a finite time step using a 
measured wind field. The growth and buoyancy of the puffs 
are computed from simultaneous specifications of atmospheric 
turbulence intensity and stability and from buoyant heat 
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flux at the source. in inversion cap through which pollu- 
tants cannot pass and the source height where pollutants are 
released are variable and can easily be adjusted. Grid dis- 
tances within the model may vary from meters to kilometers 
and time durations from seconds to hours are possible. 

This puff model has the capability of monitoring a maxi- 
mum of twenty-five sources of puffs and its grid may contain 
up to 100 puffs. A puff source can be located anywhere on 
the grid and have a unique release rate, starx and stop of 
release time, and heat production. When the center of a 
puff moves outside the boundaries of the grid (either hori- 
zontally or vertically) , that particular puff is dropped 
from memory. In this way the model doss not store irrele- 
vant puff information, thus keeping computer memory require- 
ments to a minimum. 

A variable to control the amount of reflection/absorp- 
tion of the pollutant by the surface is easily adjusted in 
the puff model. Such a parameter is of great value both in 
actual dispersion problems and also for gaining understand- 
ings of the plume/ surf ace relationship. 

The model calculates the concentration at each grid 
point by summing the contributions from surrounding puffs 
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for each advection step. The grid concentrations can 
be allowed to accumulate or simply be updated with the lat- 
est instantaneous value. A minimum grid concentration of 
interest can be set to reduce computer run time by dropping 
concentrations too small to be of interest. 

The output of the model contains periodic results of 
puff locations and concentrations as well as initial input 
verification. The time interval for the periodic results is 
adjusted by the input data. This recurrent lineprinter out- 
put contains: 

• X-Y plane plots showing the position 
of the sources and of puffs inside the 
grid, 

• X-Z plane plots of puff positions 
for evaluating plume rise for each 
vertical level of interest, and 

• a table listing of the grid point 
concentrations for each level. 

A computer drawn contour chart of the magnitudes of the pol- 
lutant concentrations is also available. 
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When considering distance between gridpoints (delta 
X,Y,Z), only spatial resolution and computer resources need 
be considered. Calculated concentration accuracy is not 
related to the grid-point separation. To insure that no 
essential information on individual puffs is "hidden” 
between grid points, the grid separation should be adjusted 
dependent upon the size of the puffs at the downwind dis- 
tance of interest. Other specific model configuration con- 
siderations are described in the following sections. They 
are also discussed in more detail in the model behavior 
chapter. 

B. HIND FIELD 

Once a puff is released, it is advected based upon wind 
data measurements at a single point only, normally the 
release point. This limits the validity of the model to 
situations where the wind field and turbulence can be 
assumed to be horizonally homogeneous throughout the grid. 
It is therefore important to insure that the data obtained 
from such a single point measurement is representative of 
the wind structure for the whole area of interest. 

The wind data are normally obtained in the form of a 
horizonal velocity time series. A vector sequence is formed 
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by averaging over a convenient interval. These data are 
read into the model after being segregated into turbulence 
classes as discussed in the next section. 

C. TOHBOLEMCE INTENSITY AND DIPFOSION 

The growth/diffusion of a puff depends upon the turbu- 
lence intensity. To account for this growth, the puff model 
applies the theory for relative diffusion suggested by Smith 
and Hay (1961). 

The turbulence intensity is defined to be the standard 
deviation (sigma) of the wind direction (in radians) 
squared. The sigma values are collected for the same short 
time periods as the wind speed measurements used to advect 
the puffs. Therefore, the intensity of the turbulence which 
governs the relative diffusion of the puffs, can be adjusted 
along with the the advecting wind speed after each time 
step, if conditions warrent. 

A very low value of turbulence intensity (as 0.0002) 
represents a small standard deviation (0.9), normally a sta- 
ble atmosphere and a weak puff dispersion/diffusion. As the 
atmosphere becomes more unstable, the turbulence intensity 
increases along with an increase of sigma values and plume 
dispersion/diffusion. While these characteristics are 
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rasentative of turbulence over land, they can be applied to 
over water cases in a broad sense. 

D. PLUHE RISE 

In the vertical direction, puff-risa can be accounted 
for by Briggs (1970) plume rise theory. In this case buoy- 
ancy is assumed to be conserved (adiabatic motion) , and 
pressure forces, molecular viscosity and local density 
changes are considered small and are neglected. The rate at 
which a puff rises as it is adv acted dowawind is a function 
of the buoyancy flux, wind speed, puff distance traveled and 
stability of the atmosphere. Plume rise is considered sep- 
erately for each individual puff. 

E. REFLECTION 

The interaction of the pollutant with the surface is 
adjustable and can be easily changed in the input data. 
Total reflection or absorption ora fraction between the two 
can be used. 

F. LIMIT OF MIXING DEPTH 

The effect of an atmospheric lid (inversion) can be 
applied in the model to limit the vertical movement of the 
pollutant. The model does not permit the plume to rise 
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above this cap. When a maximum mixing level is in effect, 
it acts to totally reflect the pollutants in the same manner 
as total reflection at the surface. This mixing cap would 
act as an inversion when the puff would be expected to grow 
much more readily in the hori zonal than in the vertical 
direction. 
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III. DMA CO LLEC TION 



An intensive field tracer study was psrformed during the 
fall of 1980 and winter of 1981 in the Santa Barbara Channel 
area of the California coast. The work was supported by the 
Bureau of Land Management (BLM) and performed by the Envior- 
mental Physics Group of the Naval Postgraduate School (NPS) 
and AeroVironment Inc (AV) , Los Angeles, CA. This study was 
designed to help validate and/or modify Gaussian dispersion 
models for coastal use and to build a data base for future 
model development. Air pollution models in current use have 
not been adequately validated for the ovec-water regime. 

In the experiment, a tracer gas (SFq) was released sev- 
eral miles offshore from the NPS Research Vessel (R/7) Aca- 
nia. Ambient gas concentrations as low as 10 parts per 
trillion (PPT) were determined by an array of land based 
sensors, from a small boat and at various levels by an air- 
craft equipped with a continuous SPg analyser. A chart of 
the experimental area and locations of the various platforms 
is shown in Fig. 2. 
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Figure 2. Expef imental Area shoyfing Locations o 
Acania, Aircraft Track and Nuraarical 



The aircraft flew through the plume at various eleva- 
tions offshore and overland. The plume transect tracks per- 
tinent to this study were made parallel tD the coast approx- 
imately one-half mile offshore. The airborne sampling 
consisted of instaneous concentrations (PPT) at selected 
points at different levels over a period of six hours. The 
observations, recorded at locations 10-70 at altitudes of 61 
and 91 m above MSL on January 29, 1980 are shown in Pig. 3. 
Average concentrations over the noted time period are shown 
at the bottom of each altitude block. 

The following marine meteorological parameters were 
measured onboard the R/7 Acania while anchored approximately 
7. 4 km offshore: 



•relative wind speed 
•wind speed fluctuation 
•sea surface temperature 
•sky cloud cover 
•relative wind direction 



•air temperature 
•dew point 
•ship roll 
♦ship location 



•inversion height (acoustic sounder) 

•vertical temperature and humidity profilies 
(shipboard radiosonde launch every 12 hours) 
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TRANSECT 
BEG END 

17:20:50-17:21:14 



17:33:35-17:34:04 



17:20:50-17:59:42 



DATE: 9/29/80 TIME: 1700-1800 PDT 

ALTITUDE: 300 FEET MSL 1/2 MILE OFFSHORE 




DISTANCE (GRID COORDINATES) 



DATE: 9/29/80 TIME: 1700-1800 PDT 

TRANSECT 
BEG END 

17:25:43-17:26:18 
17:37:00-17:37:41 
17:25:43-17:37:41 

30 40 

DISTANCE (GRID COORDINATES) 



ALTITUDE: 200 FEET MSL 1/2 MILE OFFSHORE 

1342 4967 4967 1753 707 595 184 

707 3024 5304 3921 520 558 109 

0 2129 7249 6053 3848 1419 372 166 

-75 335 5081 7922 2465 634 634 335 

671 3548 6108 3903 2277 1007 276 93 

316 1679 5192 5821 260 1511 371 317 167 

■ ‘ ' ■ ■ ■ ■ I ‘ 



Figure 3 



Aircraft observed Plume C 
Grid Coordinates at 6 1 an 
September 1980. 



oncent rations 
d 91 m above M 
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The tracer gas was released at a fixed rate through the 
exhaust of one of the ships motor generator sets. The gen- 
erator was run at a constant speed resulting in a constant 
stack temperature and flow rate. 

The above collected meteorological information enabled 
the data grid to be established and the atmospheric wind 
field, turbulence intensity, stability and buoyancy flux to 
be derived for inputs into the puff model. The puff model 
prediction based on this actual data formed the basis for 
the model performance evaluations carried out here. From 
this basis, different input variables were adjusted to note 
the effect on the advected concentrations — both in relation 
to each other and to the airborne measurements. 
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17. aODBL B EHA7 I0B 



The Puff Model was run on the NPS IBM S/370 Model 3033 
AP computer with two goals in mind; (1) familiarization with 
model performance under actual conditions and (2) a compari- 
son of model predictions with observed data. The two goals 
are interrelated in the sense that atmospheric data col- 
lected in the aforementioned tracer study were used to form 
an initial prediction of the pLume dispersion, and varia- 
tions of that data were used to evaluate the limits of the 
model. The data used as. input to the model represented the 
marine atmospheric conditions as determined from R/7 Acania 
meteorological data at the time of the experiment. 

Proper grid spacing was arrived at by considering puff 
spread, mean wind direction and the geographical area of 
interest. With the initial prediction in hand, data input 
variables of the model were adjusted and their effects 
(changes in prediction) noted. All model predictions were 
compared with the aircraft observed data. 

A 7. 4 X 4. 3 km downwind area of interest was initially 
gridded into an 17 X 10 array. Distances between horizonal 
and vertical grid points were approximately 435 m (Fig. 2). 
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since many of the aircraft observation times centered 
around 1730 hours (all times are Pacific Daylight Time), 
model puff releases were initiated at 1630. The 30-minute 
wind speed averages obtained from data taken onboard the 
ship between 1630 and 1730 were 4.7 and 4.8 ms“‘. The first 
advected puffs would be expected to arrive at the back edge 
of the grid slightly before 1700, and by 1730 a steady con- 
sistent plume would be passing through the area of the air- 
craft track. (The model showed, in fact, puffs leaving the 
back edge of the grid slightly before 30 minutes after puff 
release) . 

The average wind direction was recorded on the ship 
every 15 seconds. The standard deviations of the wind 
direction (sigma) were computed as approximate one minute 
averages. These, in turn, were averaged over 30 minutes to 
correspond with the 30-minute wind speed averages. The 
sigma values during the time of interest were 1.0 and 0.9 
resulting in the very small turbulence intensity values of 
.0003 and .0002. 

A delta Z value of 33 m was used to observe plume con- 
centrations at the altitudes of 0, 33, 66 and 99 m above the 
surface. These levels were chosed for comparison with the 
aircraft transect altitudes of 6 1 and 91 m. 
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Fine scale vertical temperature and humidity plots were 
drawn based upon radiosonde soundings taken onboard the 
ship. The sounding taken at 1735 PDT (Figure 4) shows a 
shallow unstable layer near the surface topped by an inver- 
sion extending to near 400 m. A 80 m depth of the mixing 
layer was subjectively established. The potential tempera- 
ture gradient computed by the formula 

|i = II t .0098 z (1) 

was found to be 1.0 deg K/100 m. 

Basic data to determine source strength and heat emis- 
sion from the ships stack were taken from Schacher, et al 
(1981). As previously mentioned, SFg gas was released 
through the ships motor generator exhaust at a constant 
rate. The stack temperature was 250 deg F, the flow rate 
was 7.13 ^ 103 m^s”*. The SFg release rate was 47.91 lb 
hr-i. The top of the ships stack — coasidered to be the 
source elevation — was 4m. The source strength was con- 
verted to 6.04 gm s“‘ for input into the model. 

Heat emission (H) in KW was determined by the formula 
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H 



( 2 ) 



= AT * ^ * Cp * Flow 



where 

delta T = temperature (stack - air) dag K 
E = dry air gas constant 
= 2. 87^? 10* erg/g deg K 
p = 1Q3 P(mb) = dyne/cm2 
Cp = specific heat of dry air 
= . 24 cal/g deg K 
Flow = 16.39 » 7.13 ^ 10 3 cm 3/g 

The heat emission was thus computed as 15.07 KM. 

Initially the model grid was established after noting 
the area of maximum airborne sampled concentrations (between 
points 24 and 43 of Fig. 3) and the location of the ship. 
It became obvious during early model runs that, with the 
actual wind direction input, the model predicted plume was 
being advected south of the grid towards point 60 on the 
aircraft track. Obviously, the steering wind, as measured 
onboard ship, was not constant all the way to the shore. A 
northward turning of the plums was indeed detected several 
times during the experiments by the aircraft. To compensate 
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for this effect, the source of the plume release was moved 
in the model three grid spaces ( 1305 m) to the north so that 
the maximum predicted plume concentrations would pass 
through the areas of the maximum airborne measured concen- 
trations. No corrections were made to the model predicted 
plume concentrations because of this adjustment. However, 
one could reason that the predicted concentration values 
would be higher in comparison with measured values since the 
coastal turbulence and wind shif t — which would tend to dif- 
fuse the plume--were not considered. 

Initial model runs with the small turbulence intensity 
classifications of .0003 and .0002 failed to show plume con- 
centrations greater than 1 ^ 10” ‘2 in the grid at any level 
other than at the source. Apparently, the grid spacing was 
too large and the narrow plume was advecting between the 
grid points. In an effort to locate the plume, a combina- 
tion of model runs were performed varying the turbulence 
intensity and grid spacing as shown in Table I. 

In this table, a mixing level cap of 80 m was in effect 
for the model predictions. No concentrations above that 
level were allowed in the computations. As previously men- 
tioned, the model mixing level cap totally reflects all 
pollutants bade downward. 
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with a grid spacing of 435 m and a turbulence intensity 
of .05, a plume concentration covering four grid spaces at 
the east end of the grid was produced. Predicted plume con- 
centrations slowly decreased as the turbulence intensity was 
reduced to .02 (atmospheric stability increased). At an 
intensity of .01, the concentration dropped by about five 
orders of magnitude. Normally, one would expect increased 
concentrations with increased atmospheric stability. Per- 
haps the result noted here is due to the plume shrinking 
away from a grid point (and becoming more concentrated 
between the recorded grids) with the increase in stability. 

To increase the grid resolution, the grid spacing was 
reduced by half to 217.5 m and again to 108.75 m. With each 
reduction the grid was reduced by half in the "Y” direction 
and doubled in the **X” direction thus keeping distances 
between grid spaces equal in all directions. This of course 
greatly increases the computational requirements. If only 
plume predictions along the back edge are needed (as in 
Table I) , the downwind grid distance may be held constant at 
435 m while the horizonal crosswind resolution is increased. 
In this way many unnecessary computations are not made. 
However, the increased horizonal crosswind resolution is 
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computed over the entire downwind grid, which in this case, 
is not necessary. A more satisfying solution to this prob- 
lem is to install the capability of using a variable resolu- 
tion grid with the model so that downwind areas of 
particular interest can be covered with a dense grid while 
other areas of not so much interest can be sparsely grided. 

In order for the advection of the plume to remain on the 
array when increasing the horizonal crosswind resolution and 
decreasing the area exposed on the grid, the plume source 
was adjusted along the western boundary of the grid. The 
relationship of the vertical grid points to changes of the 
source location is shown in Table II. The plume source for 
each grid resolution is noted with an arrowhead. Grid 
points that are aligned vertically in the table have identi- 
cal locations and should have the same predicted plume con- 
centrations. As mentioned earlier, an increase of grid 
resolution does not affect the predicted concentration. 
Notice that for the same grid points in Table I, the pre- 
dicted concentrations with a turbulence intensity of .01 
remain constant with changes in grid spacing — only the grid 
resolution was changed. 
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TABLE II 



The Relationship of the Y Axis along the Western Grid Edge 
to changes of Grid Distance between 435 and 108.75 m. 



DELTA 
y (m) 

435 

217.5 

108.75 



Grid Points Along the Y Axis 
4 5 6 7 8 9 

0123456789 
0123456789 

k 



Prom Table I it is obvious that the 435 m grid spacing is 
too large and that the higher resolution does indeed "see” 
concentrations that would otherwise be missed. 

The problem of a increasingly narrow distance covered on 
the grid as resolution is increased can sometimes be at 
least partly corrected by reversing the K and Y coordinates 
and adjusting or rotating the advecting wind direction. 
This can easily be done with the use of the "TURN” model 
input parameter. This procedure sometimes becomes necessary 
since one of the grid directions is limited by the width of 
the output printer paper to less than or equal 10 grid 
units. 

The model input variables, meteorological and source 
values were adjusted to note their effect on plume concen- 
trations. A deeper understanding of how the model worics and 
how the atmosphere affects dispersion can also be gained by 
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such adjustments. A turbulence intensit 7 class of .05 was 
used, except when studying intensity itself, because it had 
previously demonstrated a good downwind grid coverage of the 
plume. 

In order to note the effect of the maximum mixing level 
on the plume concentrations, several model predictions were 
run, varying only the height to which the plume was allowed 
to rise. Exact grid point reproductions were not possible 
since the model only allows the height of the mixing level 
to be an integer multiple of delta Z. The vertical grid 
spacing is therefore not equal. However, the antipicated 
trend of increased concentrations as tha mixing level is 
lowered is evident from Table III. 

The reflection/absorption of the plume at the surface is 
controlled by the model variable '‘REFLE3". Tests of the 
extremes of total absorption (0.0) and total reflection 
(1.0) were performed. The results showed a 50 percent 
reduction in plume concentratio ns at the west end of the 
grid with total absorption compared to total reflection in 
otherwise identical model runs. 

The model has a self-imposed limitation of 100 puffs 
from all sources on the grid. The model will terminate if 
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TABLE III 



Plume Concentrations between 5ui7face and 99 m under 
Different Maximum Mixing Levels. 



Max Mixing 



Level 


(5) 


(6) 


None 


.44E-8 


.55E-5 




.43E-8 


.54E-5 




.40E-8 


.50E-5 




.36E-8 


.45E-5 


80 ra 


.57E-8 


.72E-5 




.59E-8 


.74E-5 




.58E-8 


.72E-5 




0 


0 


30 m 


.64E-8 


.81E-5 




.65E-8 


. 81E-5 




0 


0 



(7) 


(8) 




20E-5 


.22E-9 


Sfc 


19E-5 


.21E-9 


33 m 


18E-5 


.20E-9 


66 m 


16E-5 


.18E-9 


99 m 


26E-5 


.28E-9 


Sfc 


27E-5 


.29E-9 


40 m 


26E-5 


.29E-9 


80 m 


0 


0 


> 80 m 


29E-5 


.32E-9 


Sfc 


29E-5 


.32E-9 


30 m 


0 


0 


> 30 ra 



this number is exceeded. A balance must be made between the 
rate at which puffs are released from the source (TAO) and 
the time it takes the puffs to be advected across the grid. 
A release rate of one puff every 40 seconds was predomi- 
nantly used during this study. 

The turbulence intensity variable was varied to include 
conditions that are more unstable. As atmospheric instabil- 
ity increases, the plume would be expected to expand 
whereas, with stable conditions, the plume should remain 
narrow and highly concentrated. 
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with the use of NPS contouring routines and the 
subroutine "DRAW”, a visual comparison of the plume disposi- 
tion and concentrations is available. Since the computed 
plume concentration varies over many orders of magnitude, 
the concentration values were converted to integer numbers 
by multiplying by 1. ^ 10^3 and then taking the logarithm. 
These logarithms are then smoothed. Thus, a contour plot 
representing order of magnitude concentrations was be pro- 
duced. As with the model variables "HAPTIH" and "KPLANS" 
which control the frequency and vertical levels of printer 
plots, "DRAW” can be called to contour concentrations at any 
time period and for any level required. 

Plume concentration distributions foe turbulence inten- 
sities of .05, .10 and .25, all other variables constant, 
are shown in Figs. 5, 6 and 7. As expected, the plume 
becomes wider and the concentration decreases as the turbu- 
lence intensity/diffusion increases and the atmosphere 
becomes less stable. 

In Figs. 5-7, the plume source was located at grid point 
(0,6). The wide plume in the lower part of the plot is not 
real but is a function of the smoothing routine spreading 
the early puffs more than would be expected. Since the 
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Figure 5. Orders of MagnitU(Je of Plume Concentration with 
Turbulence Intensity equal to .05. 
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Figure 6 
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Figure 7. Same as Figure 5 except Turbulence Intensity 
equal to .25. ^ 
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smoothing routine would tend to smooth strong concentrations 
near the source, the smoothing should be eliminated if the 
primary interest is near the source. Actually one would 
expect the puffs to behave as in Fig. 3, from Mikkelsen 
(1979), showing the relationship between the puff size and 
concentration, the rate of puff release (TACJ) and the 
advecting wind speed 0. 





Figure 8. Helationship between Puff Size, Concentration, 

Puff Release Rate (rAO) and Advecting Wind Speed 
0 (Mikkelsen, 1979) . 
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puffs would have to travel the distance Xmin before they 
expand to a size where they effectively overlap and form a 
solid plume. From Figs. 5-7, one can see that plume concen- 
trations have increased with distance and that a Xmin has 
been reached in the middle to upper part of the grid. It is 
at this point that the puff model would be expected to accu- 
rately predict plume concentrations. If the area of inter- 
est is before the present Xmin, the release rate of puffs 
would need to be increased so that the successive puffs 
would overlap sooner. Also noted that as the turbulence 
intensity increases, the area of maximum concentration of 
the plume expands while the central concentration decreases. 
This agrees with conservation of mass theorey. 

To appreciate the relative importance of the source 
strength and buoyant heat flux, these variables (discussed 
in Chapter IV) were doubled separately and together and the 
concentrations compared to the concentrations from the 
actual conditions. Little or no changes in concentration 
were noted when the buoyant heat flux was doubled and source 
strength remained the same. However, when the source 
strength was doubled and heat flux held constant, the grid 
concentrations doubled as expected. Thus, under existing 
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conditions, the source strength was critical to the pre- 
dicted plume concentrations while the buoyant heat flux, 
within the range tested, was not relevant. The vertical 
printer plots did show an initial puff rise soon after 
release due to the initial heat release but as the puff rose 
and expanded, it soon reached the ambient temperature and 
leveled off. The buoyant heat flux would probably be more 
important when dealing with a smaller scale grid or greater 
heat release. 
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V 



DATi COMPARISON 



No attempt was made to compare actual puff model concen- 
tration predictions at exact grid points to aircraft obser- 
vations for the following reasons: 

• The aircraft locations were approximations — the exact 
locations were not known. Large differences in predicted 
concentrations are seen with small grid separations as 
evidenced in Table II. 

• As noted in Fig. 2 , the aircraft observations were taken 
over a period of time at different levels — while the puff 
model produced multilevel instantaneous predictions. 

• As mentioned earlier, the actual wind was not constant 
between the ship observation site and the oppsite side 
of the grid near shore. Since the model advects the 
puffs based upon ship observed wind, the behavior of 
actual plume would be different from predicted. 

• Calabration procedures for the SF^ continuous analyzer 
mounted onboard the aircraft were not available for 
instantaneous concentrations greater than 1010 PPT. 



45 



Therefore, a question of actual levels of SFs 
concentration in the higher ranges exists. 

The puff model predicted concentrations are expressed in 
while the aircraft observations are shown as the volume 
of SFgper unit volume of air in PPT. A conversion between 
the predicted and observed concentrations was obtained by 
computing the partial pressure and molecular weight of SF© 
at standard pressure and temperature. A conversion of 
g/ffl3 = (.63*510-11) <5 observed concentration (PPT) 
was thus found. 

Generally, the aircraft sampled concentrations (Pig. 2) 
show values between 100 and 8000 PPT. Converting these 
observed concentrations to predicted concentration units 
gives values between .63*^10-« and .50510“'^ g/m^. These 
observed concentrations are much smaller than the values 
shown in Table II. Perhaps this difference could be 
explained by the fact that the puff model advected the plume 
toward the coast in the same direction under the same very 
stable conditions as observed on the ship. Any considera- 
tion of increased turbulence a nd wind shifts near shore 
would be expected to reduce the actual plume concentrations 
toward the observed concentration levels. 
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Increasing the turbulence intensity to 0.25 and keeping 
all ether variables constant, the concentration values would 
decrease to the order of magnitude of 10-* — closer to the 
observed concentrations. (Phis would require the wind 
direction standard deviation to increase from 1 to 28). 
However the increased instability would rause the plume to 
spread over a much greater area (Fig. 7) than observed by 
the aircraft. 
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VI. CONCLUSIONS MD HECOMHEND ATIONS 



The puff model has been demonstrated to be a versatile 
working dispersion model. Different combinations of input 
variables showed the expected reasonable results. The dif- 
ferences between model predicted and aircraft observed plume 
concentrations do not seem to be the fault of the model but 
mainly that of the highly variable meteorological conditions 
found along a coast. 

Probably the most obvious conclusion reached from this 
study is that predicting the behavior of a plume moving over 
a marine environment onto a coastal region has significant 
problems. In all probability, atmospheric boundary layer 
conditions 7.4 km offshore can be very different from those 
observed in the more turbulent coastal region. The single 
point meteorological measurement at the source should not be 
expected to adequately represent plume characteristics as it 
nears a meteorologically variable coastline. Additional 
observations (primarily wind spaed and direction), or other 
means of predicting the coastal meteorological conditions, 
would have to be incorporated into the puff model to ade- 
quately handle this problem. 



48 



The advantage of incorporating variable grid spacing 
within the puff model and the obvious beaifits have already 
been discussed. 

Presently, the mixing cap of the puff model is required 
to be located at an integer multiple of delta Z. More flex- 
iability in this parameter to include any level, regardless 
of delta Z, would be benifical. 

along with the puff locations shown on the lineprinter 
output, a maximum concentration level of each puff would be 
helpful. 

In future experiments, several aircraft tracks should be 
made further out from the coast in an attempt to avoid the 
turbulent coastal region. Observations thus obtained in a 
noncoastal environment would help to verify the model pre- 
dictions without the coastal influence. 
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APPENDIX A 



MAJOR SECTIONS OF THE POFF 33DEL 

The Ris0 Puff Model has been described by (Mikkelsen, 
1979). The code also is well documented with comment state- 
ments. With that information and the outline to be provided 

/ 

in this and the following appendices, the computational and 
input/output procedures will be obvious. 

The program and input data are st<^ei on cards for the 
sake of permanency. For efficient operational execution, 
the program and input data cards are read on a disk within 
the computer. The model can then be run at will without 
reference to the original data cards. Minor changes can 
easily be made directly on the disk both to the model and/or 
data before each execution. 

The model can be separated into the following main sec- 
tions: 

A. Input Data 

B. Initial 

c. Calculating 

D. Output 

E. Error Diagnostics 

F. Subroutines 
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These will be described separately in the following sec- 
tions. 

a. INPOT DATA SECTION 

The input data includes the variables shown in Table IV. 

TABLE IV 

Input Data Variables for the Puff Dispersion Model. 

Mind History Potential Temperature Gradient 

Turbulence Intensity Buoyant Heat Flux 

Grid Dimensions Minimum Cnncentration of Interest 

Mixing Depth Reflection at Ground Level 

Source Locations, Start/Stop Time, Strength, Heat Emission 
Number of Seconds between Advection Steps 
Number of Seconds between Pr intouts/PLots 
Number of Seconds between Puff Releases 



The wind field and stability class for the current rime step 
are read at the start of the calculation section. 

The variables listed above are printed as a input data 
check and a permanent record to accompany the actual output. 
In most cases the print command can be o/erridden by YES/NO 
options. 

B. INITIAL SECTION 

Based upon the input data from section (A) , the initial 
section specifies and initializes parameters to be used in 
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the calculating section and is passed only once during exe- 
cution of the model. The grid and some counters are ini- 
tialized. Constants relating to reflectance, mixing depth 
and stability as well as those controlling the size of some 
of the loops within the model are established. Parameters 
such as number of puff releases per second, number of advec- 
tion steps per second and number of advection steps per puff 
release are determined. 

C. CALCOLATION SECTION 

Osing current wind and stability class data read at the 
start of the calculation section, the model advects the puff 
centers and calculates the growth rate and plume rise of the 
puffs. It removes the puffs that have left the grid (hori- 
zonally and/or vertically) . The predicted concentration is 
computed at the grid points to include pollutants from all 
nearby puffs. 

D. OUTPUT SECTION 

For time intervals designated by the input data, printer 
plots of the X-Y and Y-Z grid are produced. A maximum mix- 
ing level is marked on the Y-2 grid if in effect. 
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These plots include the source location and a trace of the 
plume from the release time to the maptime. Also printed at 
this interval is a X-Y table of grid concentrations for each 
vertical level of interest. These concentrations can be 
either accumulated or actual concentrations at the plot 
time. 

Added to the puff model is a versatic plotter routine to 
smooth and contour the grid magnitude concentrations of the 
above tables. 

E. ERROR DIAGNOSTIC SECTION 

If the model is directed by the input data beyond the 
limits of the design of the program, the program is termi- 
nated by way of the error diagnostic section. It prints 
comments relating to commonly made input errors enabling the 
user to isolate problems. 

F. SOBRODTINES 

The subroutine '’Sigris” calculates the puff size in the 
horizonal and vertical directions. It also estimates plume 
rise associated with pollutent buoyancy. 

The subroutines "Ispace" and "R space” are used in the 
framework of the printer plots. 
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The subroutine "Draw” converts the plume concentrations 



to a logarithmic values, smoothes and then contours 
using NPS inhouse contour subroutines. The values are 
verted to their logarithm values so that the problem of 
touring over many orders of magnitude is simplified. 



them 

con- 

con- 
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APPENDIX B 



PUFF MODEL FLOW CHART 



INPUT 

SECTION 



READ - ciata input variables (primary, source, 

intensity, current stability) ; averaging 
■ timestep (tav) for wind data 



INITIAL SECTION [ Initialize constants 



READ & STORE all wind speed and direction 
for current intensity class 



CALCULATION 



next wind vector 



SECTION 



yes 



I RECALL f rcm storage one wind vector for next tav 

next advection step 



Release puff, execute advection step, 
calculate concentration 







E ‘advection^ 
steps >-no 
.GE. tay,< 



A 



E ^vection 
^ steps 
.BQ. MAPTIM 



.yes. 




OUTPUT 

SECTION 

Puff positions, 
grid concentra- 
tions, contour 



no 



fEnd ol 


datal 







I Termination! 
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APPENDIX C 



POFF MODEL CONTRACTIONS 

CHEMIN — Minimum grid concentration of interest 
DELX ,DELY, DELZ — Distance in meters between grid points 
DOSE — Allows the concentration matrix to accumulate 
DTDZ — Potential temperature gradient (K/M| (.GE. 0) 

HEAT — Individual source heat emission (KW) 

ICOLS — Number of columns in grid (.LE. 10) 

INST — Instantaneous concentration matrix 

ITIME — Start time 

JROWS — Humber of rows in grid 

KPLANS — Number of vertical levels in grid (includes surface) 
MAPTIM — Number of seconds between printer plots 
NRELSE — Number of seconds ro stop of release 
NRMDLT — Number of sources (.LE. 25) 

NTADV — Integer number of seconds between advection steps 
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REFLEC — Reflection at ground level (0. none; 1.0 total) 



SOORNR — Number to identify source 

SOORST — Strength for individual source (gm/s) 

STOPRL — Individual source stop time (s) 

STRTRL — Individual source start time (s) 

TAO — Integer number of seconds between puff releases 
TORN — Angle of rotation of wind dirsction 
XSOORCE--X coordinate of source in grid units 
YSOOBCE--Y coordinate of source in grid units 
ZH--Limited mixing depth (m) 
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